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Abstract 

It is known that classical string dynamics on pure AdS$ x 5 5 is integrable and plays 
an important role in solvability. This is a deep and central issue in holography. Here we 
investigate similar classical integrability for a more realistic confining background and 
provide a negative answer. The dynamics of a class of simple string configurations on 
AdS soliton background can be mapped to the dynamics of a set of non-linearly coupled 
oscillators. In a suitable limit of small fluctuations we discuss a quasi-periodic analytic 
solution of the system. Numerics indicates chaotic behavior as the fluctuations are not 
small. Integrability implies the existence of a regular foliation of the phase space by 
invariant manifolds. Our numerics shows how this nice foliation structure is eventually 
lost due to chaotic motion. We also verify a positive Lyapunov index for chaotic orbits. 
Our dynamics is roughly similar to other known nonintegrable coupled oscillator systems 
like Henon-Heiles equations. 
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1 Introduction 



AdS /CFT duality pQ was a major step towards the goal of recasting large- iV QCD as a string 
theory. One of the interesting aspects of AdS /CFT is integrability (see the review [2]). Integra- 
bility has allowed us to obtain many classical solutions of the theory that would otherwise have 
been impossible to find [3], 0]. One may approach integrability from two sides corresponding 
to the two extreme values of the 't Hooft coupling. On the supergravity side (which is a good 
description as A — > oo), integrability of the classical sigma model on AdS$ x S 5 was established 
for bosonic sector in [5] and fully completed with the inclusion of fermions in [6]. It has been 
shown that classical string motion in AdS*, x S 5 has an infinite number of conserved charges^] 
This is the closest to solvability that we can get currently. It is also to be noted that the study 
of integrability in non-linear sigma models has a long history [S] • The other approach to 
integrability is perturbative in the weakly coupled gauge theory [2J. We will restrict ourselves 
to the former approach. 

An important open question is whether integrability can be extended to more QCD-like 
theories]^] The original form of AdS /CFT duality was for conformally invariant J\f = 4 SYM 
theory but this can be deformed in various ways to produce string duals to confining gauge 
theories with less or no supersymmetry. The construction of [SI E] does not readily generalize 
to these less symmetric backgrounds. One prime example of a confining background is the 
AdS soliton [9j [10]. Similar geometries have been used extensively to model various aspects 
of QCD in the context of holography [11] . Here we will look at the question of integrability of 
bosonic strings on an AdS soliton background. Although it is much more interesting to explore 
full quantum integrability, to begin with we may ask whether we can find enough conserved 
charges even at a purely classical level. The answer turns out to be negative. By choosing a 
class of simple classical string configurations, we show that the Lagrangian reduces to a set 
of coupled harmonic and anharmonic oscillators that correspond to the size fluctuation and 
the center of mass of motion of the string. The oscillators decouple in the low energy limit. 
With increasing energy the oscillators become nonlinearly coupled. Many such systems are well 
known to be chaotic and nonintegrable [T2l IT3] . It is no surprise that our system also shows 
a similar behaviour. Possibly chaotic behaviour of a test string has been argued previously in 
black hole backgrounds [H] [15] . However our problem is somewhat different as we are looking 
at a zero temperature geometry without a horizon. In a companion paper [16] non-integrability 
of string theory in AdS$ x T 1,1 is discussed. 

We start by discussing the AdS soliton background and our test string ansatzS. We 
discuss some quasi-periodic solutions for small oscillation regime. Our argument for the non- 
integrability is via numerical solution of the EOM's^H In a certain regime of parameter space 
the system shows a zigzag aperiodic motion characteristic of a chaotic system. We then look 
at the phase space. Integrability implies the existence of a regular foliation of the phase 
space by invariant manifolds, known as KAM (Kolmogorov-Arnold-Moser) tori, such that 
the Hamiltonian vector fields associated with the invariants of the foliation span the tangent 
distribution. Our numerics shows how this nice foliation structure is gradually lost as we 

4 However the commutator algebra of such charges is not fully understood. 
5 This is one of the motivations discussed in the introduction of [B]. 
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increase the energy of the system^! To be complete we also calculate Lyapunov indices for 
various parameter ranges and find large positive values in chaotic regimes^^H We discuss open 
questions and possible extensions in the conclusion^! 



2 Setup 

The AdS soliton (Ai) metric for an asymptotically AdSd+2 background is given by [9], 

ds 2 = L 2 a' ( e 2u (-dt 2 + T 2 Ju)d9 2 + dw 2 ) + ^^du 2 
V " T 2n (u) 

where T 27T (u) = 1 - I ^-eM . (1) 

At large u, T 2n (u) ~ 1 and (JIJ reduces to AdSd+2 in Poincare coordinates. However one 
of the spatial boundary coordinates 9 is compactified on a circle. The remaining boundary 
coordinates Wi and t remain non-compact. The dual boundary theory may be thought of as a 
Scherk-Schwarz compactification on the 9 cycle. The 9 cycle shrinks to zero at a finite value of 
u, smoothly cutting off the IR region of AdS. This cutoff dynamically generates a mass scale 
in the theory, very much like in real QCD. The resulting theory is confining and has a mass 
gap. 

Here we will work with0 d = 4 and make a coordinate transformation u = Uq + ax 2 with 
u = log(2/5) and a = 5/4, such that T 27T (u ) = and for small x the x-9 part of the metric 
looks flat, ds 2 ~ dx 2 + x 2 d9 2 . In these coordinates the metric is 




dg 2 = e 2u 0+ 2a X ^_ dt 2 + j.^2 + rf 2) + 

T(x) 

where T(x) = 1 - e~ 5ax2 . (2) 



2.1 Classical string in AiS-soliton 

We start with the Polyakov action: 

Sp = ~2^ J d ' rdff V = 7'f b G llv d a X' t d b X v (3) 

where X M are the coordinates of the string, G^ u is the spacetime metric of the fixed background, 
jab is the worldsheet metric, the indices a, b represent the coordinates on the worldsheet of the 
string which we denote as (r, a). We work in the conformal gauge jab = Vab and use the 
following embedding for a closed string (partially motivated by [14"]): 

t = t(r), 9 = 6(t), x = x(t), 

W\ = R(t) cos (4>(cr)) , w 2 = R(t) sin (4>(cr)) with 0(a) = aa . (4) 
6 The analysis for any other d > 3 proceeds along the same lines and almost identical results can be obtained. 
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The string is at located at a certain value of u and is wrapped around a pair of w-directions 
as a circle of radius R. It is allowed to move along the potential in u direction and change its 
radius R. Here a G Z is the winding number of the string. The test string Lagrangian takes 
the form: 

" {-i 2 + T(x)6 2 + w 2 -w' 2 } + ^x 2 (5) 



C oc -e 
5 



2ax* 



T(x] 



'2a.v 



j-t 2 + T(x)6 2 + R 2 - R 2 a 2 j + 



2^2 



2a z x 



-x 



(6) 



where dot and prime denote derivatives w.r.t r and a respectively. The coordinates t and 6 are 
ignorable and the corresponding momenta are constants of motion. The test string Lagrangian 
differs from a test particle Lagrangian because of the potential term in R(r). The coordinate 
R would be ignorable without a potential term. In general it can be easily argued that for a 
generic motion of a test particle in an AdS soliton background, all the coordinates other than 
x are ignorable and the equations of motion can be reduced to a Lagrangian dynamics in one 
variable x. This implies integrability. 

Here the conserved momenta conjugate to t and 6 are, 

Pt = --e 2ax2 i = -E 



Pe 



e 2ax2 T(x)9 = k. 



(7) 



The conjugate momenta corresponding to the other coordinates are: 

4 

Aa 2 x 2 



PR 



: e 2ax R 



Px 



T(x) 



■x . 



With these we can construct the Hamiltonian density: 



U = 



8 



-E 2 + 



k' 2 



T(x) 



+ Pr 



-2ax 2 _|_ T(x)p 2 x 16 j}2^2 a 2ax 2 



Hamilton's equations of motion give: 
R 



\PR^ 



PR 



X 



--Ra 2 e 2ax2 
5 

T(x)p x 



4a 2 x 2 



(8) 

(9) 

(10) 

(11) 
(12) 



Px = 



Aax 



m*)pi 

5a 2 x 3 



+ 



E 2 



Pi 
5a 2 x 2 



T(x) 



P\ 



e -2ax- + l A R 2 a 2 e 2ax 2 

25 



k 2 e 
T(x) 



-2ax l 



d x T{x) 



(13) 
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We also have the constraint equations: 

{d T X»d T X v + d a X»d a X v ) = , (14) 
G^drX^X" = . (15) 

The first equation takes the form H = Q and the second equation is automatically satisfied 
for our embedding. 



3 Dynamics of the system 

At k = 0, an exact solution to the EOM's is a fluctuating string at the tip of the geometry, 
given by 

x{r) = (16) 
R(t) = Asin(r + <p). (17) 

where A, <j) are integration constants. No such solution with constant x(t) exists for k ^ 0. 
However one may construct approximate quasi-periodic solutions for small R(t),pr(t). It 



should be noted that with R,pr = the zero energy condition Eqn.(14) becomes similar to 
the condition for a massless particle and the string escapes from AdS following a null geodesic. 
For small nonzero values of Ro,pn, the motion in the x-direction will have a long time period. 
However the fluctuations in the radius will have a frequency proportional to the winding 
number which is of 0(1). This is a perfect setup to do a two scale analysis. In the equation 
for pn we may replace R{t) 2 by a time average value. With this approximation, motion in 
the x-direction becomes an anharmonic problem in one variable which is solvable in principle. 
The motion is also periodic [Fig |l(a)j . On the other hand to solve for R(r) we treat x(r) as a 
slowly varying field. In this approximation the solution for R(t) is given by 

R(t) « exp(-a x(r) 2 )A sin(r + 0). (18) 



Hence R(t) is quasi-periodic [Fig, 1(a)] . We have verified that in the small R regime, the 



semi-analytic solution matches quite well with our numerics. 

Once we start moving away from the small R limit the the above two scale analysis breaks 
down and the nonlinear coupling between two oscillators gradually becomes important. In 
short the coupling between oscillators tends to increase as we increase the energy of the string. 
Due to the nonlinearity, the fluctuations in the x- and .R-coordinates influence each other and 
the motions in both coordinates become aperiodic. Eventually the system becomes completely 
chaotic [Fig, l(c)| . The power spectrum changes from peaked to noisy as chaos sets in [FigjTj . 



As we discuss in the next subsection, the pattern follows general expectations from the KAM 
theorem. 



7 The Hamiltonian constraint could be tuned to a nonzero value by adding a momentum in a decoupled 
compact direction. For example if the space is Ai x S 5 then giving a non-zero angular momentum in an S 5 
direction would do the job. However we choose to confine the motion within M. here. 
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E=0.22, k=0.2, x(0)=l.l 




Frequency 

(a) Periodic motion showing R(t) ,x(t). (b) Power spectrum of x(r) for periodic motion. 
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(c) Chaotic motion of string showing x(t). (d) Power spectrum of x(t) for chaotic motion. 

Figure 1: Numerical simulation of the motion of the string and the corresponding power spectra 
for small and large values of E. The initial momenta Px(0),p#(0) have been set to zero. For a 
small value of E = 0.22, we see a (quasi-)periodicity in the oscillations. The power spectrum 
shows peaks at discrete harmonic frequencies. However for a larger value of E = 3.0, the 
motion is no longer periodic. We only show x(r) but R{t) is similar. The power spectrum is 
white. 

3.1 Poincare sections and the KAM theorem 

An integrable system has the same number of conserved quantities as degrees of freedom. A 
convenient way to understand these conserved charges is by looking at the phase space. Let 
us assume that we have a system with N position variables with conjugate momenta p^. 
The phase space is 2iV-dimensional. Integrability means that there are N conserved charges 
Qi = fi{P) <?) which are constants of motion. One of them is the energy. These charges define a 
iV-dimensional surface in the phase space which is a topological torus (KAM torus). The 2N- 
dimensional phase space is nicely foliated by these iV-dimensional tori. In terms of action-angle 
variables (Jj, 9i) these tori just become surfaces of constant action. With each torus there are 
N associated frequencies uJi(Ii), which are the frequencies of motion in each of the action-angle 
directions. 
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It is interesting to study what happens to these tori when an integrable Hamiltonian is 
perturbed by a small nonintegrable piece. The KAM theorem states that most tori survive, 
but suffer a small deformation [T2"l [T5] . However the resonant tori which have rational ratios 
of frequencies, i.e. rriiUi = with m G Q, get destroyed and motion on them become chaotic. 
For small values of the nonintegrable perturbations, these chaotic regions span a very small 
portion of the phase space and are not readily noticeable in a numerical study. As the strength 
of the nonintegrable interaction increases, more tori gradually get destroyed. A nicely foliated 
picture of the phase space is no longer applicable and the trajectories freely explore the entire 
phase space with energy as the only constraint. In such cases the motion is completely chaotic. 

To numerically investigate this gradual disappearance of foliation we look at the Poincare 
sections. For our system, the phase space has four variables x, R,p x ,pR. If we fix the energy 
we are in a three dimensional subspace. Now if we start with some initial condition and time- 
evolve, the motion is confined to a two dimensional torus for the integrable case. This 2d torus 
intersects the R = hyperplane at a circle. Taking repeated snapshots of the system as it 
crosses R = and plotting the value of (x,p x ), we can reconstruct this circle. Furthermore 
varying the initial conditions (in particular we set R(0) = 0, p x {0) = 0, vary x(0) and determine 
Pr{0) from the energy constraint), we can expect to get the foliation structure typical of an 
integrable system. 

Indeed we see that for smaller value of energies, a distinct foliation structure exists in the 
phase space [Fig. 3(a)| . However as we increase the energy some tori get gradually dissolved 

[Figs 



mm 

3(c)P(d) 



The tori which are destroyed sometimes get broken down into smaller tori 
Eventually the tori disappear and become a collection of scattered points 
known as cantori. However the breadths of these cantori are restricted by the undissolved tori 



and other dynamical elements. Usually they do not span the whole phase space [Figs. 3(c)- 



3(f)] . For sufficiently large values of energy there are no well defined tori. In this case phase 



space trajectories are all jumbled up and trajectories with very different initial conditions come 
arbitrary close to each other [Fig, 3(h)] , The mechanism is very similar to what happens in 



well known nonintegrable systems like Henon-Heiles models 
3.2 Lyapunov exponent 

One of the trademark signatures of chaos is the sensitive dependence on initial conditions, 
which means that for any point X in the phase space, there is (at least) one point arbitrarily 
close to X that diverges from X. The separation between the two is also a function of the initial 
location and has the form AX(Xo, r). The Lyapunov exponent is a quantity that characterizes 
the rate of separation of such infinitesimally close trajectories. Formally it is defined as, 

A= lim lim *i A*(*o,r) 

r^ooAXo^or AX(X ,0) v ; 

In practice we use an algorithm by Sprott [17], which calculates A over short intervals and 
then takes a time average. We should expect to observe that, as time r is increased, A settles 
down to oscillate around a given value. For trajectories belonging to the KAM tori, A is zero, 
whereas it is expected to be non-zero for a chaotic orbit. We verify such expectations for our 
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case. We calculate A with various initial conditions and parameters. For apparently chaotic 
orbits we observe a nicely convergent positive A [Fig|2] . 



E=0.22, k=0.2, x(0)=l.l E=3,k=0.2, x(0)=l.l 




(a) Lyapunov index for periodic motion. (b) Lyapunov index for chaotic motion. 



Figure 2: Lyapunov indices for the same values of parameters as in Fig.([l|. For E = 0.22, 
the Lyapunov exponent falls off to zero. For E = 3.0, the Lyapunov exponent converges to a 
positive value of about 0.38. 



4 Conclusion 

In this work, we argue using numerical techniques that classical string motion in the AdS 
soliton background is nonintegrable. This certainly restricts the solvability of such theories. 
Also our results give a perspective on how much of classical integrability may be extended 
to various holographic backgrounds, especially those with less symmetry than AdS§. Non- 
integrability is possibly quite generic and might be demonstrated by studying time evolution 
of simple string configurations. In particular, the basic construction of our work seems to be 
extendible to other confining backgrounds [18J. It would be nice to explore these directions. 
Another interesting extension would be to include world sheet fermions. It is also to be kept in 
mind that the AdS soliton is not an exact string background and possibly has a' corrections. 
However these effects are unlikely to change the main result of the current work. 

One big question is the implication of our result for the full quantum spectrum. This is 
in turn connected with the glueball spectrum of the dual theory. The low-lying string modes 
will possibly be decoupled from the center of mass motion and will be more like the flat space 
counterpart. However the higher modes will be affected by the nonlinearity. Many exact results 
are known for the quantum spectrum of a chaotic theory. It would be interesting to explore 
how these results apply in a mini-superspace quantization of our system. 

It would also be interesting understand more on the gauge theory side [19]. The dual 
gauge theory is the = 4 SYM theory with one compact direction with aperiodic boundary 
condition for fermions. This breaks supersymmetry and the low energy dynamics of the theory 
is confining. Here, a simple change in boundary condition is changing the integrability of the 
theory. It is also to be noted that the full SYM theory with 4> corrections is almost surely 
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nonintegrable. Any apparent integrability would then be a property of the large-iV saddle 
points. 
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(a) E = 0.25 



(b) E = 0.30 
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(g) E = 1.0 



(h) E = 3.0 



Figure 3: Poincare sections demonstrate breaking of the KAM tori en route to chaos. Each 
colour represents a different initial condition. For smaller values of E the sections of the KAM 
tori are intact curves, except for the resonant ones. The tori near the resonant ones start 
breaking as E is increased. For very large values of E all the colours get mixed - this indicates 
that all the tori get broken and they fill the entire phase space. 
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